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Based on exact diagonalization and variational cluster approximation calculations we study the 
relationship between charge transfer models and the corresponding single band Hubbard models. 
We present an explanation for the waterfall phenomenon observed in angle resolved photoemission 
spectroscopy (ARPES) on cuprate superconductors. The phenomenon is due to the destructive 
interference between the phases of the 02p orbitals belonging to a given Zhang-Rice singlet and 
the Bloch phases of the photohole which occurs in certain regions of k-space. It therefore may be 
viewed as a direct experimental visualisation of the Zhang-Rice construction of an effective single 
band model for the Cu02 plane. 

PACS numbers: 74.72.-h,79.60.-i,71.10.Fd 



I. INTRODUCTION 

The 'waterfall phenomenon' observed in angle resolved 
photoemission spectroscopy on cuprate superconductors 
has attracted some attention. This phrase summarizes 
the following phenomenology [H-fTijj : for photon energies 
around 20 eV, where photoholes are created predomi- 
nantly in 02p statesflj], a 'quasiparticle band' can be 
seen dispersing away from the Fermi energy as T is ap- 
proached but then - along the (1, 1) direction roughly at 
(f i \ ) - rapidly looses intensity and cannot be resolved 
anymore. Instead there is a 'band' of weak intensity 
which seems to drop almost vertically in k space towards 
higher binding energy, thus creating the impression of a 
kink in the quasiparticle band. Up to binding energies of 
« 1.0 eV there is a 'black region' around the T-point with 
no spectral weight at all. Finally, at binding energies of 
around 1.0 eV the 'waterfall' seems to merge with one or 
several bands of high intensity, which often correspond 
very well to bands predicted by LDA calculations. 
In experiments with photon energies ~ lOOeV (where an 
appreciable fraction of the photoholes is created in Cu3d 
states [HI) and when the spectra are taken in higher Bril- 
louin zones fA Hoj|. however, the low energy quasiparticle 
band can be resolved all the way to T and shows no indi- 
cation of a kink. This band therefore undoubtedly exists 
and has no kink so that the only possible explanation for 
the 'black region' and the apparent kink seen at low pho- 
ton energies are matrix element effects. This conclusion 
has in fact been reached by Inosov et al. on the basis of 
their experimental data(7l fioj. A very similar conclusion 
was also reached by Zhang et aZ.[13J who showed that 
the kink in the quasiparticle band appears only in the 
second derivative of the momentum distribution curves 
but is absent in the second derivative of the energy distri- 
bution curves which shows a smooth quasiparticle band 
instead. They concluded that the kink is not an intrinsic 
band dispersion. 

Some support for this point of view comes from the fact 
that the waterfall phenomenon is observed over the whole 
doping range, from the antiferromagnetic insulator to the 



Fermi-liquid-like overdoped compounds, so that it is ob- 
viously unrelated to any special features of the electronic 
structure. Moreover, the spectrum of a hole in an anti- 
ferromagnetic insulator is one of the very few reasonably 
well-understood problems in connection with cuprate su- 
perconductors. For this special problem good agree- 
ment between experiment [la. Il6j . various approximate 
calculation s 1 1 7H24j . and exact diagonalization of small 
clusters [25l - l28j has been found and no theory for hole 
motion in an antifcrromagnet predicts a kink in the quasi- 
particle band. Here in particular the work of Leung and 
Gooding[27j should be mentioned who studied the spec- 
tral function for a single hole in a 32-site cluster by ex- 
act diagonalization and found a well-defined and smooth 
quasiparticle band in excellent agreement with the re- 
sults of the self-consistent Born approximation [3, l45j - 
Moreover, since plasmons with energies below the charge 
transfer gap are not expected in the insulating com- 
pound and since the coupling to spin excitations is ob- 
viously described very well by the self-consistent Born 
approximation, a hypothetic Bosonic mode can be al- 
most certainly ruled out as an explanation of the water- 
fall phenomenon in the undoped compounds and, due to 
the near-independence of the phenomenon on the doping 
level, also for the entire doping range. 
One type of matrix-element effect which partly explains 
the nonobservation of the quasiparticle band at T has 
been discussed by Ronning et al. [1| . At the T-point in the 
first Brillouin zone the photoelectrons are emitted exactly 
perpendicular to the Cu02 plane which we take as the 
(x — ?/)-plane of the coordinate system. In this situation - 
and if we neglect small deviations from this symmetry in 
the actual crystal structure - the experimental setup has 
Ci v symmetry. The expression for the photocurrent [39| 
involves the dipole matrix element (f\A • p|i), where A 
is the vector potential of the incoming light, p is the mo- 
mentum operator, |i) the initial and |/) the final state. 
The state |/) differs from \i) by 

a) the presence of an electron in the so-called LEED 
state, which far from the surface evolves into a plane 
wave oc e tkz and hence transforms according to the iden- 



tical representation 

b) by the presence of a hole with momentum (0,0) in 
a Zhang- Rice singlet. Since the Zhang-Rice singlet has 
a Cn3d x 2_ y 2 orbital as its 'nucleus' [40j it has the same 
symmetry. 

It follows that |/) and \i) have the same parity under 
reflection in the x — y plane so that the dipole matrix 
element is zero if A is in the CuC>2 plane. On the other 
hand |/) and \i) aquire a relative minus-sign under a ro- 
tation by \ around the z-axis so that the dipole matrix 
element is zero as well if A is perpendicular to the Cu02 
plane. This means that it is impossible to observe a state 
with the character of a Zhang-Rice singlet at T in nor- 
mal emission. As Ronning et al. pointed out, however, 
this argument cannot explain the nonobservation of the 
quasiparticle band at T in higher Brillouin zones, where 
the photoelectrons are no longer emitted in the direction 
perpendicular to the surface. 

Moreover, a similar effect has been observed in ARPES 
studies of the compound SrCuC>2 which contains CuO 
chains. With a photon energy of 22.4 eV and polar- 
ization parallel to the chains - which generates holes in 
cr-bonding 02p orbital - there is no intensity at k = o(29j. 
If the photon energy is increased to 100 eV, however, the 
spinon-band around k = can indeed be resolved (30j. 
This behaviour is similar to the waterfall effect but in this 
compound the Cu02 plaquettes are perpendicular to the 
surface of the crystal. We conclude that there must be 
a second mechanism for the extinction of spectral weight 
around T in both, the one and two dimensional systems. 
An explanation of the waterfall phenomenon has been 
given by Basak et a/.[3l|. These authors first showed 
that the waterfall phenomenon cannot be reproduced by 
a calculation of ARPES spectra from LDA band struc- 
tures and eigenfunctions alone, a procedure which has 
otherwise been found to be highly successful in describ- 
ing ARPES spectra of cuprate superconductors [321 [33|. 
Instead, these authors obtained good agreement with ex- 
periment by additionally introducing a self-energy which 
describes the coupling to a Bosonic mode. The basic 
mechanism for the waterfall and the variation of the 
ARPES spectra with photon energy then is bilayer split- 
ting which is enhanced by the coupling to the Bosonic 
mode, in particular the vertical part of the waterfall turns 
out to be the strongly renormalized bonding combination 
of the two single-layer wave functions. This model repro- 
duces the strong changes of ARPES spectra with photon 
energy hv in the range 60 eV — > 80 eV as observed by 
Inosov et aZ.0, E3| quite well. On the other hand, In- 
osov et al. gave a different explanation for this variation, 
namely the rapid variation of the Cu3d photoemission in- 
tensity at the Cu 3p — > 3c? absorption threshold at 75 eV. 
Since the strong variation of photoemission spectra at the 
3p — s- 3d threshold is well established for 3d transition 
metal oxides [lij this appears a plausible explanation. 
There have been attempts to reproduce a kink in the 
band structure within the framework of a single-band 
Hubbard or t-J model(35Tj37|. The high- intensity bands 
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FIG. 1: (Color online) Possible realization of the model {TJ. 

observed near T at binding energies of ~ 1.0 eV thereby 
are identified with high ener gy features observed previ- 
ously in exact diagonalization |25l - l27| or Quantum Monte- 
Carlo 38] studies for such models. 

In the present manuscript we take the point of view that 
the waterfall phenomenon is a pure matrix-element ef- 
fect, as pointed out by Inosov et al. 0, E3] and Zhang et 
al.^^. Thereby the crucial point is the nature of the 
low-energy hole states as Zhang-Rice singlets. Since the 
phases of the 02p hole 'within' a Zhang- Rice singlet cor- 
respond to momentum (-7T, 7r)[40] there is perfect destruc- 
tive interference with the phases of a p-like photoholc 
with momentum (2wr, 2mn) with n and m integer. The 
quasiparticle band at T therefore can be observed only at 
photon energies where the cross section for hole creation 
in d-orbitals is large - because the destructive interference 
occurs only for 02p-photoholes - and in higher Brillouin 
zones where the argument by Ronning et aL[l[ does not 
apply. 

In section II we will discuss exact diagonalization results 
for a 1-dimensional charge transfer model. In section III 
we discuss the spectra of a 2-dimensional charge transfer 
model by the variational cluster approximation (VCA). 
In section IV we discuss the experimental relevance of 
the binding energy of the Zhang-Rice singlet and section 
V gives summary and conclusions. 



II. ONE DIMENSIONAL MODEL 

We study a 1-dimensional (ID) charge transfer model 
by exact diagonalization. We choose a ID model because 
we need at least a two-band model and the largest clus- 
ter of a two-band model we can study in 2D contains 4 
Cu-ions, so that we have virtually no k-resolution. As 
will be seen below, however, a very simple ID model is 
sufficient to reproduce the key features of the waterfall 
phenomenon. To be more precise, the Hamiltonian reads 

H = tY,[4,APi +ht ,-Pi- h „) + H.c: 

i,cr 

+A ^ n d,i + U ^ n dM n d,i,l- (1) 

i i 

where d\ a (Pj a ) creates a hole in a d-like orbital at site i 
(p-like orbital at site j). We choose t as the unit of energy 
and use the values A = — 4 and U = 8 throughout. The 
crucial feature of the model is that - as in the case for 
the Cu02 plane- the d — p hybridization integral has an 
alternating sign. A schematic representation of the model 
is shown in Figure [T] Figure [H shows the single particle 
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FIG. 2: (Color online) Single particle spectral function for 
the model ([TJ for the half-filled system (8 holes in 8 unit 
cells). The Figure combines spectra calculated with different 
boundary conditions to give an impression of larger systems. 
The right part of the spectra corresponds to hole creation 
(i.e. photoemission) the left part to hole annihilation (inverse 
photoemission). The upper figure shows the spectra for U — 
0, the lower one the spectra for (7 = 8. 



spectral function for the noninteracting case U = and 
for the strongly correlated case U — 8. These are denned 
as 
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3 G dd (k,cj + irf) 



(2) 



where Gdd is the d-like diagonal element of the 2 x 
2 single-particle Green's function (and analogously for 
A p (k,u))). The Lorentzian broadening r/ = 0.1. In 
the noninteracting case there are two bands, one with 
predominant p-character and one with predominant d- 
character. Since the matrix element of the p — d hy- 
bridization is oc sin(-|) there is no p — d hybridization for 
k = and the states have pure p- or d-character. Also, 
the energy of the p-like peak agrees exactly with that of 
the p-orbital, i.e. 0. 

Surprisingly the p-like band can be identified also in the 
strongly correlated case, particularly so near k = 0. The 
reason is that the p-like Bloch state with k = does not 
mix with the d-like Bloch state due to parity so that a 
hole created in this state is unaffected by the Coulomb 
repulsion on the d-sites. Accordingly, the energy of this 
peak still is 0. This state is analogous to the T eV-peaks' 
in the real cuprate materials [52]. Near k — the mixing 
still is small so that U is effectively only a weak pertur- 
bation. 

On the other hand the upper band of predominant d- 
character disappears completely and is replaced by a 
spectrum which is very similar to that of a ID single- 
band Hubbard model. This can be seen in Figure [3] which 
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FIG. 3: (Color online) Top: Closeup of the hole addition part 
of the spectrum in Figure [5] Bottom: Photoemission spec- 
trum for a half-filled single-band Hubbard chain with 8 sites 
and U'/t' = 10.8 and t' = 0.335 t. The single-band spec- 
trum has been turned 'upside down' to be compatible with 
the hole picture and shifted by 0.87t. The figure combines 
spectra calculated with periodic and antiperiodic boundary 
conditions. 



shows a close-up of the low energy region of the photoe- 
mission (i.e. hole addition) spectrum and compares this 
to the photoemisssion spectrum of a single-band Hub- 
bard model. It turns out that for the above values of t, 
U and A a very good match can be obtained by choosing 
the parameters of the single-band Hubbard model to be 
U'/t' = 10.8 and t' = 0.335 t. There is a rather obvious 
one-to-one correspondence between the peaks, the differ- 
ent 'holon bands' characteristic for the spectra of finite 
clusters of the one-dimensional Hubbard or t-J model 43 1 
can be identified in both spectra. The main difference 
occurs at energies E between —0.5 and 0. This is most 
likely the consequence of mixing and level repulsion be- 
tween the single-band Hubbard-like bands and the free- 
electron-like p band. Moreover one can see a splitting of 
some of the peaks in the two-band model. 
The full interacting spectrum therefore can be modelled 
well by superposing the lower noninteracting p band and 
a single band Hubbard spectrum. This indicates that 
the Zhang-Rice construction [40] for reduction of the low- 
energy sector of the two-band model to a single band 
model also works well for this ID model. On the other 
hand, for higher energies the 2-band model has additional 
states which correspond to the non-bonding combination 
of p-orbitals. 
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FIG. 4: (Color online) Closeup of the hole addition part of 
the spectrum in Figure [2] The left panel shows the p-like 
spectrum, the right panel the d-like spectrum. 



Figure |4] shows the photoemission spectrum of the two- 
band model split into its p- and ci-component. This 
demonstrates that the key features @, [HI of the water- 
fall phenomenon can be seen already in this simple ID 
model: Near T there is no p-like intensity in the single- 
band Hubbard-like states, instead all p-like intensity re- 
sides in the free-electron-like band (this is the huge peak 
at E = 0). The p-like intensity alone therefore shows 
the black region around k — and a 'high-energy-kink' 
whereby the vertical part of the waterfall is due to the 
incoherent continuum of the single-particle spectrum of 
the single-band Hubbard model. In the d-like spectrum 
on the other hand the quasiparticle band can be fol- 
lowed right up to the T-point. Already this simple model 
therefore reproduces the key features of the waterfall phe- 
nomenon and the photon energy dependence of the spec- 
tra: the waterfall occurs in the first Brillouin zone for all 
photon energies and in higher Brillouin zones for photon 
energies where the cross section for Cu3d is small. 
To understand the extinction of p-like weight around T 
we repeat the Zhang-Rice construction [|(| and consider 
a single plaquette - which in the ID model consists of the 
(i-orbital and its two nearest neighbor p-orbitals - at the 
d-site Hi. For one hole, the ground state reads 

l*o,J = («i + «i 4 )|0) 

- >l^-^J- (3) 

The coefficients {u\,V\) are the GS eigenvector of the 
matrix 



The singlet ground state - i.e. the analogue of a Zhang- 
Rice singlet in the ID model - of two holes is 

l*^> = ( «a pV-,i + ^ (dW~,i+P-A ) 

+w 2 44 )|o) 

and the coefficients (u 2 , v 2l w 2 ) are the GS eigenvector of 
the matrix 

/ 2t \ 
H 2h = \ 2t A 2t . (5) 
\Q 2t 2A + U J 

In the corresponding single-band model, the state l^o ff ) 
corresponds to the site i being occupied by a spin-cr elec- 
tron, whereas \^q) corresponds to a hole at site i. The 
matrix element of the electron annihilation operator Ck,-|- 
between theses states is e lkRi . On the other hand, the 
matrix elements of the operators pk.-j- and dk.t are 

k 1 
m p {k) = e lkRi V2 isin(-) (mu 2 + —^viv 2 ) 

2 v2 

m d (k) = e ik - Rl (-^=u lV2 +V!W 2 ). (6) 
V 2 

The crucial term here is the factor of i sin(|) which arises 
from the overlap of the p-like Bloch state with momen- 
tum k and the bonding combination pj _ a . In addi- 
tion, we have to take into account a shift in energy of 
e = Eq' 1 — Eq. This corresponds to the binding energy 
of the ZRS and has no counterpart in the single-band 
model. In simplest terms, we would therefore expect 
that the photoemission part of single particle spectral 
functions of the two-band model can be obtained from 
that of the single band model by 

A p , d {k,uj) = \m p , d (k)\ 2 A(k,uj + e). (7) 

In this expression several simplifications have been made: 
the ZRS is assumed to extend only over one plaquette, 
which is probably not correct. This implies that pro- 
cesses where a ZRS in a plaquette around site i is gen- 
erated by actually creating a hole in a neighboring unit 
cell are neglected. Moreover the problem of the overlap of 
Zhang-Rice singlets in neighboring cells [13] is not taken 
into account either. 

The energy shift e which is necessary to match the two 
spectra in Figure [3] is found to be 0.87 t - the estimate 
obtained from the eigenvalues Eq H and E 2h is 1.04 t, i.e. 
reasonably close. Figure [5] shows the spectra obtained 
from equation ([7]) compared to the actual spectra of the 
two-band model. The numerical values of the prefactors 
are uiu 2 + -^y\v 2 = 0.70 and -^uiv 2 + vitu 2 — 0.50. 
While the agreement is not really perfect, the qualitative 
trends are reproduced well, particularly so near k = 
where the hybridization with the p-like band is weak. 
The main differences occur for larger fc-values and may 
also be due to the fact that the ZRS is not restricted to 
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FIG. 5: (Color online) Closeup of the photoemission spectrum 
in Figure [2] (Top) compared to the spectra of a single-band 
Hubbard model corrected according to ([7]) (Bottom). The 
energy shift e is 0.84 t rather than 1.04 t as would be obtained 
from the single-plaquette calculation. 



one plaquette and also the hybridization with the p-like 
band which is absent in the single-band Hubbard model. 
On the other hand, given the simplicity of the proce- 
dure for converting the single-band spectra into two-band 
spectra this not so bad and qualitatively explains the 
extinction of p-like intensity around V at least qualita- 
tively: this is due to the factor of sin(|) in the matrix 
element m p , which describes the destructive interference 
between the phase factors of the two p-orbitals in the 
bonding combination p[_ a and in the electron operator 

p^ k . This is in turn the consequence of the oscillating 
sign of the hopping integral in (JTJ) so that the same mech- 
anism should also be effective in the 2D CuCVplane. 
The exact diagonalization results then can be summa- 
rized as follows: The oscillating sign of the d — p hy- 
bridization induces an oscillating sign also in the bond- 
ing combination of p-orbitals around a given d-site. The 
phase of the p-orbitals in the bonding combination there- 
fore 'locally' corresponds to a momentum of it. It follows 
that around k — 2mr there is destructive interference 
between the phases of the bonding combination and the 
phases of the photohole, and it is not possible to couple 
to the ZRS by hole creation in p-orbitals. The ZRS-part 
of the spectrum thus becomes extinct in the p-like spec- 
trum, whereas no such extinction occurs in the rf-likc 
spectrum. To conclude this section we note that none 
of the above considerations is limited to half-filling. Fig- 
ure [HI compares the spectra for the 2-band model and the 



FIG. 6: (Color online) Low energy photoemission spectrum 
of the hole-doped two-band model compared to the spectrum 
of the corresponding single-band model. In this figure, the 
chemical potential is the zero if energy. 



single-band model in the hole-doped case, that means at 
a hole density of 1.25. Thereby all parameters are the 
same as in Figure 03 Again, a very good correspondence 
between the two models exists and again the extinction 
of p-like weight around k — can be clearly seen. 



III. 2 DIMENSIONAL MODEL 

We now apply the picture gained from the analysis of 
the one-dimensional two-band model to obtain approxi- 
mate ARPES spectra for the 2D CuCVplane. Since it is 
not possible to compute the spectra of a 2D three-band 
model larger than 2x2 unit cells by exact diagonaliza- 
tion we switch to the variational cluster approximation 
(VC A) [H, [Ho|, EH to compute at least approximate spec- 
tra. The VCA uses the fact 51] that the grand canoni- 
cal potential of an interacting Fermi system can be ex- 
pressed as a functional of the self-energy S(w) and is 
stationary with respect to variations of at the ex- 

act self-energy. The VCA then uses finite clusters - the 
so-called reference system - to numerically generate 'trial 
self energies' for an infinite system. This is described in 
detail in the literature [i^, [50, [H[ and has turned out to 
be a very successful method to discuss the single-particle 
spectral functions of correlated electron systems. We use 
this method to calculate spectra for the three-band Hub- 
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FIG. 7: (Color online) Top: Sum of p-like and d-like spectrum 
for the three-band model ((SJ with U = 0. Center: d-like 
spectrum for (7 = 8 as obtained by the VCA. The spectrum 
has been multiplied by 2.5. Bottom: p-like weight for U = 
8. The spectra are computed for half-filling, that means 1 
hole/unit cell, the black line denotes the respective chemical 
potential. 
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a=x j£L a ,a 

-A ^2 n d ,i + U d d ^ n d,i,t n dM 
ieL d i&L d 

+ U pp ^ n p,j,t n p,j,l- 

jeL a 



(8) 



Ld denotes the s.c. lattice of Cu3d sites, L x and L y the 
s.c. lattices of p x and p y sites. Moreover 



K = - (p' 

2 V x 



■ x — P ■ , » —P.« + P^ ■ , v 

I— §,<r r x,i+^,a C y,i—%,(T y,i+Q,cr 



(9) 



is the bonding combination of p-orbitals around the d- 
orbital at site i and 
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for a — x denotes the bonding combination of p y orbitals 
around a given p x orbital and vice versa for a — y. a de- 
notes the unit vector in a-direction. The model is again 
formulated in hole language, i.e. eft creates a hole in 
orbital i and t p d, t pp > 0. We choose t p d as the unit of en- 
ergy, the other parameters are Udd = 8, U pp = 3, A = 3, 
t pp = 0.5. We study this model by the VCA, using a 
cluster with 2x2 unit cells (i.e. a square shaped CU4O8 
cluster) with 4 holes (corresponding to 'half filling') as the 
reference system for creating trial self-energies. Some re- 
sults obtained in this way for this model have previously 
been published by Arrigoni et aL53j. Figure [7] shows the 
total spectral weight for the case case U = as well as 
the p-like and d-like spectral weight for the interacting 
case, U — 8. The p-like spectral weight now is defined as 
the sum of the two p-like diagonal elements of the total 
3x3 spectral weight matrix. 

In the noninteracting case, U — 0, there are three bands. 
Switching on the Coulomb repulsion has a very similar 
effect as for the ID model: The two upper (in hole lan- 
guage) bands with predominant p character remain essen- 
tially unchanged, whereas the partly filled band of pre- 
dominant d-character is split into an upper and a lower 
Hubbard band. Comparing the p and d spectra in the 
energy range — > 2t p d the same extinction of p-like spec- 
tral weight around (0, 0) can be seen i.e. exactly the same 
behaviour as in the ID model. 

By analogy with the ID case we assume that the reason 
for the extinction of p-like weight again is the matrix el- 
ement between a plane wave of p holes and the bonding 
combination ©. We again consider a single-plaquette 
problem. For later reference we include a Coulomb re- 
pulsion Upd between p and d-holes on neighboring sites 
which is set equal to zero for the time being. The Hamil- 
ton matrices for the single and two hole plaquette prob- 
lems are: 
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2tpp 2tpd 



2t 



pd 
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21, 




% 2y/2t 



pd 



A 



- 2/ 

^Lpp 

2y/2t 



U, 



pd 



pd 





2\/2tpd 
2A + Udd 

The p — p Coulomb repulsion oc U pp is treated in mean- 
field theory. The ground state eigenvectors of these ma- 
trices are again denoted by (ui,vi) and (1L2, V2, u^)- 
To discuss the p-like photoemission spectrum we need the 
matrix element between the bonding combination ([§]) and 
a p-like Bloch wave. We write this as 
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FIG. 8: (Color online) Right part: d-like spectrum of the 
three-band model (top) compared to the spectrum of the sin- 
gle band Hubbard model multiplied by 0.5 (bottom). 
Left part: p-like spectral weight of the three-band model mul- 
tiplied by 2 (top) compared to the spectrum of the single band 
Hubbard model multiplied by 0.5 /(k) in (I14[) (bottom). 



where m a (a € x, y) are matrix elements for hole creation 
in a p a orbital at the origin, which in a real experiment 
depend e.g. on the photon polarization and wave vector 
of the photoclcctrons. In the present calculation A p (k, u) 
is obtained by calculating the spectra for creating holes 
in p x and p y orbitals separately and adding them - ac- 
cordingly, the k-dependent correction factor in ([7]) should 
be replaced bv[54| 



77lp(k) = {UXU2 + —j^lxVi) 2 /(k) 
2 / k x 



/(k) 



sin 



(14) 



To check this we have calclated the photoemission spec- 
trum of a the single-band Hubbard model by the VCA 
using again a 2 x 2 cluster as reference system. We use the 
2x2 cluster as reference system in order to make the spec- 
tra of the two models as comparable as possible. By using 
the same reference system gaps in the VCA bands which 
originate from the supercell structure which is inavoid- 
ably introduced by the VCA should be more or less iden- 
tical in both spectra. The parameter values of the single 
band model are again chosen to obtain a good match with 
the dispersion of the three-band model. Thereby in ad- 
dition to the nearest neighbor hopping integral t' also 
hopping integrals between second-nearest neighbors 
and t' 2 between third nearest neighbots were included. 
A good match was obtained by using t'/t p d = 0.38, 
t[/t' = -0.145, t' 2 /t' = 0.118 and U'/t' = 10.3. Figure! 



compares the low energy hole addition spectrum of the 
three band model with the (shifted) photoemission spec- 
trum of the single band model. The d-like spectrum of 
the three-band model should be roughly identical to that 
of the single band Hubbard model and the two spectra on 
the right part of the figure are indeed quite similar. There 
are two main differences: the weakly dispersive band at 
« 1.75t p d which can be seen in the single band model is 
absent in the spectrum of the three-band model. This 
may be a consequence of hybridization with the two ad- 
ditional p-likc bands in the three-band model. Moreover, 
the spectrum for the three-band model has a less smooth 
dispersion with additional gaps along (0,0) — > (it,0). 
This is likely due to the lower symmetry of the three- 
band model: in the single band Hubbard model, the 2x2 
cluster with open boundary conditions is equivalent to a 
4-site chain with periodic boundary conditions. For the 
two band model a similar reduction is possible, but with 
two p orbitals in between any two successive d orbitals. 
The resulting degeneracy of the ligands in the 4-site chain 
may give rise to the additional band splitting. The left 
part of the Figure compares the p-like spectrum of the 
three-band model and the single band model spectrum 
corrected by the factor m p (k). Again, the dispersion of 
the spectral weight along comparable bands is very sim- 
ilar in the two spectra. The energy shift to allign the 
single-band and three-band spectra is e = 2.00 t p d- The 
estimate obtained from the ground state energies of the 
matrices ([11]) and (HU) is Etf 1 - Ef 1 = 2.238 t pd . The 
error of ~ 10% seems reasonable taking into account the 
various approximations made. To match the intensities 
of the single-band and three-band spectra in Figure [8] 
the single band spectra were multiplied by a factor of 
0.5. The estimates for the correction factors from the 
single-plaquette problem are {u\U2 + -^y\^i) 2 = 0.56 

and (-^uiW2 + W1W2) 2 = 0.31. 

To summarize the discussion so far: the low-energy hole 
addition spectrum of the three-band model can be ob- 
tained to reasonable approximation from an effective 
single-band model whereby the p-like spectrum needs 
to be corrected by a simple k-dependent factor which 
originates from the interference between the phases of p- 
orbitals in the bonding combinations and the phases in 
a Bloch state with momentum k. If this correction is 
done p-like and d-like spectra can be obtained to good 
approximation from he spectrum of a single band Hub- 
bard model. 

We now apply this finding to compute approximate the 
spectra for a Cu02 plane, but this time make use of the 
fact that in a single-band Hubbard model larger clusters 
can be used as a reference system. We again apply the 
VCA but thus time we use a 4 x 4 cluster with peri- 
odic boundary conditions as reference system. We used 
a cluster with 2 holes corresponding to a hole concentra- 
tion of 12.5%. It is important to use as large a cluster 
as possible for the exact diagonalization step, because 
only large clusters reproduce the incoherent continua in 
the spectral function sufficiently well and, as will be seen 




below, these incoherent continua are crucial to explain 
the experimentally observed spectra. Solution of a 4 x 4 
cluster by exact diagonalization is only possible with pe- 
riodic boundary conditions, which are not customary in 
VCA calculations. On the other hand, we do not want 
to discuss the phase diagram of the Hubbard model but 
we mainly want to obtain the spectral function so this is 
justified. In fact, the spectral function obtained by the 
simpler cluster perturbation theoryJH, [H[ is almost ex- 
actly the same as the one obtained by VCA. 
Having obtained an approximate spectrum for the single- 
band Hubbard model we again use ([7]) together with (fT4")l 
to obtain the d-like and p-like intensity for the three-band 
model. This is shown in Figure |9] for momenta along the 
(1,1) direction. 

It can be seen that the Figure reproduces the water- 
fall effect quite well. In the d-like spectrum the low- 
energy quasiparticle band can be seen together with a 
broad high intensity part at more negative binding en- 
ergy. Moreover there is appreciable incoherent weight 
around k = (0,0). These incoherent continua are 
well known from exact diagonalization [25l - l28l | and self- 
consistent Born calculations (44l l45j for the t-J model. 
Their intensity decreases with increasing distance from 
k = (0, 0). This decrease is due to the coupling of pho- 
toholes to spin and charge fluctuations [Hj]. In the p-like 
spectrum the factor m^(k) in (fl"4"|) creates the dark region 
around k = (0,0). Accordingly, the quasiparticle band 
seems to disappear at rj (^, j). Since the incoherent 
weight is reduced by the same factor m^k), it disap- 
pears at the same momentum which creates the impres- 
sion of a 'band' which has a kink at sa (f , f )• At the kink 
the spectral weight of the band drops sharply. Qualita- 
tively this is exactly what is seen in the ARPES spectra 
which show the waterfall phenomenon. The chimney-like 
appearance of the spectum is due to the fact that the d- 
like spectral weight can be approximated by a product: 
A p (k,uj) = \m p (k)\ 2 A(k,uj + e). For fixed ui the second 
factor, A(k, ui + e) decreases with |k| - see the left part of 
Figured]- whereas the first factor, |m p (fc)| 2 vanishes at 
k = (0,0) and increases with |k|. Accordingly, the p-like 
spectral weight must go through a maximum and this 
maximum corresponds to the apparent vertical part of 
the band. For a quantitative discussion it would be nec- 
essary to take into account also the interference between 
hole creation in Cu3d 2 .2_ J( 2 and 02p orbitals. This how- 
ever would necessitate to know the relative magnitude 
and phase of the respective dipole matrix elements and 
this is beyond the scope of the present paper. 



IV. DIRECT MEASUREMENT OF THE 
BINDING ENERGY OF THE ZRS 

Lastly we wish to point out that the energy shift e - i.e. 
the binding energy of the ZRS - does indeed have some 
relevance for the interpretation of experimental data and 
has in a sense been observed directly. We refer to the re- 




(0,0) 



fa*) (0,0) 



(71,71) 



FIG. 9: (Color online) Left: Spectral function of the the sin- 
gle band Hubbard model with U = 10, t = 1. 
Right: Spectral function of the the single band Hubbard 
model multiplied by /(k) (see (|14[) b 



suits of Meevasana et a/.[H) who reported an anomalous 
enhancement of the noninteracting bandwidth in Bi2201. 
These authors pointed out that the energy difference be- 
tween an assumed band bottom at F and the Fermi en- 
ergy seems to be larger than the occupied bandwidth pre- 
dicted by LDA calculations and concluded that there is 
an anomalous correlation induced band widening rather 
than the expected correlation narrowing. 
We now estimate the position of /i with respect to an 
'leV peak as an intrinsic reference energy and show that 
the Fermi energy can be obtained quite accurately from a 
single-band Hubbard or t-J model by consequent applica- 
tion of the Zhang- Rice constructionjiol . To begin with, 
Meevasana et al. observed a band at T with downward 
curvature and a binding energy of w — 1 eV relative to 
the Fermi level at T (this is the band labelled B in Figure 
1 of Ref . Q ) . The authors point out that this is an umk- 
lapp of a band at the F-point, or (71", 7r) in a simple cubic 
2D model. This is therefore probably the same state as 
shown in Fig. 3(c) of Ref. 42j, i.e. a state composed of 
02p7r orbitals which right at (71, it) has zero hybridiza- 
tion with any of the correlated Cu3d orbitals. Inspection 
of Fig. 3(c) of Ref. 42] shows that the energy of an elec- 
tron in this state at (ir, 7r) - and accordingly its umklapp 



at r - is Ei = e p + 4t 



P p where e p is the orbital energy for 



02p electrons. Since this is a single-particle-like state, 
the corresponding binding energy in the photoemission 
spectrum is Eg - E^~ l = E x . 

Next, we consider the Fermi energy. The largest part of 
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the energy shift thereby is the binding energy of the ZRS, 
e, which was discussed above. In the matrices (jlll) and 
(fT2| the energy of a hole in an 02p orbital, i.e. — e p , was 
chosen as the zero of energy so that we have to change 
e — > e + e p . It remains to add the Fermi energy /i# of the 
single-band Hubbard model itself: 

fi — e p + e + fin 

so that fjb—Ei = e+fiH — ^tpp To evaluate e we use the pa- 
rameter set given by Hybertsen et a/.[4a| in their Table I. 
The only exception is the direct oxygen-oxygen hopping. 
Here we use the value t pp = 0.37el^ which has been ex- 
tracted directly from experiment in Ref. (42j where it 
also was found to be consistent with previous estimates. 
We then obtain e = 1.82 eV. The Fermi energy of the 
single-band Hubbard model may be estimated from exact 
diagonalization results. In a 4 x 4 cluster this was found 
to be 1.6 t at U/t = 8^ and 1.788 t at U/t = 10[H]. 
Hybertsen et al. estimated U — 5.4 eV and t — 0.43 eV 
so that U/t = 12.6. We estimate (j, H « 2 t = 0.86 eV so 
that eventually fx — E\ = 1.2 eV. The value in Figure 1 
of Meevasana et al. is w 1 eV. The agreement is reason- 
able given the uncertainty about some parameters but 
it is quite obvious that taking into account the binding 
energy e is indispensable to obtain a correct estimate of 
the Fermi energy. 

The value of fi — E\ obviously depends on t pp so that 
small variations from one compound to the other may 
well explain the variations observed by Meevasana at 
al. . The Fermi energy of the Hubbard model will change 
with doping as well, but these changes are a faction of 
J m 120 meV . Accordingly, the distance between the 
free-electron-like state at T and the Fermi energy should 
always be ps 1 eV and this is indeed the case in Bi2201 
and Bi22120|. All in all one can say that the consequent 
application of the ZRS picture can explain the position 
of the Fermi energy relative to a leV peak - which forms 
a natural intrinsic reference energy - quite well. 

V. SUMMARY AND DISCUSSION 

In summary, we have investigated the relationship be- 
tween the single-particle spectra of actual charge-transfer 
models and corresponding 'effective' single band Hub- 
bard models. In the noninteracting case U = the charge 
transfer models have several bands and it was found that 
those bands with predominant ligand (i.e. p character in 
the present models) are almost unaffected by the strong 
correlations. Bands with predominant <i-character are 
split into two Hubbard bands which can be mapped quite 
well to those of an effective single-band Hubbard model. 
It turned out that the spectra of p and d electrons to 
good approximation can be obtained from that of the 
single band Hubbard model by a constant shift in energy 
- the binding energy of the Zhang-Rice singlet - and, in 
the case of the p-like spectrum, by a factor which might 
be termed the form factor of the Zhang-Rice singlet. 



These results give a natural explanation for the waterfall 
phenomenon in terms of a pure matrix element effect - as 
has previously been inferred by Inosov et al. @, E3] and 
Zhang et aZ.[l3j from an analysis of their experimental 
data. Here one has to distingish between two different 
effects: the first one - pointed out already by Ronning 
et a/.[l[ - applies only to the special situation of near- 
normal emission of photoelectrons. Here the vanishing of 
the dipole matrix element (/|A • p\i) makes it impossible 
to observe ZRS-derived states. 

The second matrix element effect is the form factor of 
the Zhang-Rice singlet mentioned above, which describes 
the interference between the phases of the p-like photo- 
hole and those of the 02p-orbitals 'within' a Zhang-Rice 
singlet - which locally correspond to a state with mo- 
mentum (-7T, 7r). This makes it impossible to couple to 
a ZRS-derived state by creating a p-like photohole with 
momenta near (2mr, 2rrnr) with integer n and m. 
These two simple rules explain under which experimen- 
tal conditions the waterfall is observed or not: in the 
hole-doped compounds and at photon energies where pre- 
dominantly 02p-like holes are produced, the quasiparti- 
cle band around Y can be observed neither in the first nor 
in any higher Brillouin zone and instead the waterfall ap- 
pears. If photon energies where Cu3d holes are generated 
are used, the quasiparticle band cannot be observed in 
the first BZ, but in higher Brillouin zones. In this case 
the waterfall is absent and no kink in the quasiparticle 
band appears. 

In the electron-doped compounds the quasiparticles cor- 
respond to extra electrons in Cu3d a; 2_ a 2 orbitals so that 
the considerations regarding the form factor of the ZRS 
do not apply. The argument regarding the vanishing of 
the dipole matrix element (/|A • p|i), however, remain 
unchanged so that the quasiparticle band around Y can- 
not be observed in the first Brillouin zone either - this 
has indeed been observed by Ikeda et aL[ll[ and Moritz 
et ai.[I3|. It should be possible, however, to observe the 
full quasiparticle band without a kink in a higher Bril- 
louin zone. 

In those cases where the 'dark region' around Y is present, 
the apparent vertical part of the waterfall corresponds 
to the incoherent continua in the single-particle spec- 
tral function of the t-J model. Since these continua are 
formed from states which also correspond to Zhang-Rice 
singlets, they become extinct near (2n7r, 2m7r) as well. 
Additional evidence comes from the fact that these band 
portions show the same dependence on photon polariza- 
tion as the quasiparticle band itself [5] . 
Finally, the high intensity bands observed at Y at binding 
energies higher than s» 1 eV have no correspondence in a 
single-band Hubbard or t-J model - otherwise they would 
not be observed in normal emission - but are precisely the 
bands of predominant 02p-character which remain unaf- 
fected by the strong correlations. They are analogous 
to the leV^-peaks observed at high-symmetry points in 
cupratesJUJ and correspond to 02p drived states which 
have little or no hybridization with the strongly corre- 
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lated d-orbitals. These states therefore are essentially 
single particle states, which immediately explains their 
much higher intensity as compared to the quasiparticle 
band. 

All in all the present theory indicates that the water- 
fall phenomenon constitutes an experimental proof of 
the Zhang-Rice construction of a single-band Hubbard 
or t-J model to describe the low energy states of the 
CuC>2 planes and in fact provides a direct visualization 
of the energy range in which the states of the real Cu0 2 - 
plane correspond to those of a single-band Hubbard or 
t-J model. It shows moreover that the incoherent con- 
tinua predicted by various calculations for the t-J or Hub- 
bard model are indeed observable in experiment in that 



they are responsible for the vertical part of the waterfalls 
themselves. 
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